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ABSTRACT 



Context. Changes in the solar surface area covered by small-scale magnetic elements are thought to cause long-term changes in 
the solar spectral irradiance, which are important for determining the impact on Earth's climate. 

Aims. To study the effect of small-scale magnetic elements on total and spectral irradiance, we derive their contrasts from 3-D 
MHD simulations of the solar atmosphere. Such calculations are necessary since measurements of small-scale flux tube contrasts 
are confined to a few wavelengths and suffer from scattered light and instrument defocus, even for space observations. 
Methods. To test the contrast calculations, we compare rms contrasts from simulations with those obtained with the broad-band 
filter imager mounted on the Solar Optical Telescope (SOT) onboard the Hinode satellite and also analyse centre-to-limb vari- 
ations (CLV) . The 3-D MHD simulations include the interaction between convection and magnetic flux tubes. They have been 
run with non-grey radiative transfer using the MURaM code. Simulations have an average vertical magnetic field of 0G, 50G, 
and 200G. Emergent intensities are calculated with the spectral synthesis code ATLAS9 and are convolved with a theoretical 
point-spread function to account for the properties of the observations' optical system. 

Results. We find reasonable agreement for simulated and observed intensity distributions in the visible continuum bands. 
Agreement is poorer for the CN and G-Bands. The analysis of the simulations exhibits a potentially more realistic centre- 
to-limb behaviour than calculations based on 1-D model atmospheres. 

Conclusions. We conclude that starting from 3-D MHD simulations represents a powerful approach to obtaining intensity con- 
trasts for a wide wavelength coverage and for different positions on the solar disk. This also paves the way for future calculations 
of facular and network contrast as a function of magnetic fluxes. 



f*^ > Key words. Sun: surface magnetism - Sun: activity 



1. Introduction 



X 



Understanding the variability of the Earth's atmosphere 
requires the ability to distinguish between anthropogenic 
and natural drivers of climate change. One important 
natural driver influencing our climate is solar irradi- 
ance (i.e. the radiative flux of the Sun, measured above 
the Earth's atmosphere) and in particular its varia- 
tion. Changes in the Sun's total irradiance influence 
the amount of energy reaching the Earth. In addition, 
changes in the solar spectral irradiance, i.e. the changes 
in the Sun's brightness at a certain wavelength, influ- 
ence _tbe^hfimistry of the Earth's upper atmosphere 
20071 ) . Finally, the h eliospheric magnetic field 
19981: ISimpsonl Il998l l shields the Earth from 
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galactic cosmic ray s which have been propos ed to in- 
crease cloud cover ( Marsh fc Svensmark 2000l) . The so- 
lar surface magnetic field, which manifests itself in, e.g., 



sunspots, faculae, and the network is able to reproduce 
solar irradiance variations on solar-cycle time scales (e.g. 
Krivova et ai1l2003t IWenzler et "al"1l2005l l2006h . For longer 



time scales, however, there is considerable uncertainty. 
While sunspots and faculae dominate the solar irradiance 
changes on timescales of days to weeks, the small-scale 
magnetic elements in the network are considered to be 
the m ain cause of long-term changes in the solar irradi- 
ance (jFoukal fc Leanlll988l ISolanki fc Fliggel 120021 1. The 
contribution from sunspots to irradiance changes is well 
understood and produces a decrease in solar irradiance. 
The contribution from faculae and mainly from smaller- 
scale magnetic elements composing the photospheric net- 
work increases the solar irradiance and is still being inves- 
tigated. In particular, the spectral response is poorly un- 
derstood. For a better understanding of irradiance varia- 
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Fig. 1. Spectrum of the visible wavelength regions at disk 
centre (calculated with ATLAS9 for the non-magnetic 
simulation) with the positions of the five Hinode filters 
indicated. 
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Fig. 2. Comparison of Hinode datasets I (Mercury transit) 
and II for the green continuum wavelength at the limb 
angle /i = 0.9. For dataset I, we plot 37 histograms where 
each includes 10 random extracts of the whole observed 
image (grey lines) and the mean histogram (solid line). We 
overplot the mean of the analogous quantity for dataset 
II (black dashed line). 

tions, it is important to study the centre-to-limb variation 
and the wavelength dependence of the emergent intensity. 

The direct measurement of the solar network contrast 
to determine the impact of these small-scale magnetic el- 
ements remains difficult, since it is (even for space in- 
struments) affected by the instrument's aperture, scat- 
tered light, instrument defocus, etc. and additionally lim- 
ited to a few wavelengths. Therefore, we calculate emer- 
gent intensities from 3-D simulations of solar magneto- 
convection and calculate the intensity variation between 
brighter and darker features to compare the simulations 
with observations. Filtergrams with a high and con- 
stant spatial resolution can be obtained with the broad- 
band filter image r (BFI) mounted on t h e Solar Optical 
Telescope (SOT, iTsuneta et all (|2008l ): llchimoto et~aT 




board the HINODE spacecraft (Ko suei et al 
a comparison of measured and simulated granulation 
contrasts, we convolved the simulations with the point- 
sprea d function (PSF) of the observation's optical sys- 
tem ( Mathew et al.l [2009) to take into account its ef- 
fects. Here, we analyse simulations with different av- 
erage magnetic field, while previ ous studies mainly in- 
cluded nonmagnetic simulations: iDanilovic et al] (|2008[ ) 
analysed Hinode/SP data and compared the inten- 
sity contrast of a Hinode/SP continuum map of the 
quie t Sun with predictions of nonmagnetic si mulations, 
and IWedemever-Bohm &: van der Voortl ( 2009h discussed 
Hinode/SOT observations in comparison to nonmagnetic 
simulations to determine the continuum intensity distri- 
butions of the solar photosphere. 

The paper is organised as follows: the observations 
are described in Sect. 2, the simulations, the intensity 
synthesis originating from them and the subsequent con- 
volution are introduced in Sect. 3. We compare the 
Hinodc/SOT/BFI observations with the simulations in 
Sect. 4 (in terms of histograms and rms contrasts) and 
also discuss the statistical aspects of this comparison. The 
centre-to-limb behaviour of the intensity contrast in our 
simulations is analysed in Sect. 5. 

2. Observations 

In this work we consider 2 different datasets from the 
Hinode Solar Optical Telescope (SOT): 

Dataset I was recorded during the Mercury transit on 
8th Nov 2006 (21:38:31-21:39:16 UT) in 5 wavelength fil- 
ters: CN bandhead ('CN', 388.3 nm), G band ('GB', 430.5 
nm), blue continuum ('BC, 450.4 nm), green continuum 
('GC', 555.0 nm) and red continuum ('RC, 668.4 nm) 
(Fig. [1]). The observed region corresponded to quiet Sun. 
The full observed Hinode image covers approximately 223 
x 112 arcsec and the pixel size is 0.05 arcsec. The helio- 
centric coordinates at the centre of the image were (0,- 
419), corresponding to a limb angle of /i = 0.9, where cos 
9 = /i, and 8 is the the angle between the line of sight and 
the line through the centre of the Sun. The dark correc- 
tion and flat fielding were done with SSW-IDL routines. 
The Mercury transit allowed the calculation of the point- 
spread function (PSF) with high confidence. Details of the 
calculation of the PSF and of thi s particular data set can 

(|2009h . 



be found in Mathew et al 



Dataset II was recorded on 8th March 2007 (18:02:32- 
20:51:46 UT) in the same wavelength filters as dataset I. 
While quiet Sun was observed also in this observation, the 
binning was different with a pixel size of 0.10 arcsec. The 
data were recorded in 12 steps from the east to the west 
limb, resulting in the availability of data at all limb angles. 
At each step, the observed image covers approximately a 
field of view of 223 x 112 arcsec, as in dataset I. For each 
image, we chose extracts such that their area corresponds 
to the desired limb angle solely. This allowed us to use 
dataset II for comparisons with calculations at limb angles 
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Fig. 3. Simulations (original and convolved) and observations in the 5 Hinode/SOT filters. Top: Original simulation 
snapshots with an input magnetic field of 50G. Black: 52% below the mean, white: 68% above the mean value. 
Middle: Same simulation snapshots convolved with an appropriate point-spread function (PSF). Bottom: Extracts 
from observed image. The black/white scale for the convolved simulations and the observations is such that black is 
33% below the mean, and white is 40% above the mean value. 
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Fig. 4. Histograms of the intensity distribution in the orig- 
inal simulation for the GC wavelength at the limb angle 
H = 0.9: 0G (black dotted), 50G (black dashed), 200G 
(black long dashes). We overplot (in grey) the correspond- 
ing histograms for the convolved simulation to illustrate 
the effect of convolving with the PSF. 



from disk centre to the limb. It should be noted that we 
are observing different areas on the Sun at different limb 
angles and that their activity levels might be different. 

Dataset I was compared with dataset II at [i — 0.9 to 
check the validity of the PSF for dataset II. An example 
is shown in Fig. [3J where we compare histograms for the 
two datasets in the green continuum (GC). The number 
of resulting subimages in the two datasets are naturally 
different: For dataset I, we segmented the observed image 
into 373 sub-images (extracts) of the desired size. Four 
extracts had to be excluded due to the Mercury transit. 
To illustrate the spread in dataset I, we overplot in Fig. [2] 
(in grey) 37 histograms where each includes 10 random 
extracts of the whole observed image. For dataset II, we 
have 151 subimages, and, for ease of comparison, we con- 
fine ourselves to plotting only the mean histogram over 
all the 151 histograms of the subimages. For dataset I, 
the subimages are scaled to the pixelsize of subimages in 
dataset II, which, for the limb angle /i = 0.9 implies 68 x 
76 pixels per subimage. The mean of dataset II appears to 
be a bit wider than that of dataset I, but overall we find 
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a good agreement. We follow the p rocedure outlined in 
Numerical Recipes (j Press et al . 1992) and calculate the x 
statistic for two binned distributions (with unequal num- 
bers of data points), where Ri is the number of events in 
bin i for the first data set, and Si the number of events in 
the same bin i for the second data set: 



x 2 = E 



R /M. 

R 



Ji V s 



(1) 



course of the 5-minute oscillations in the box. Each simu- 
lation was run for about 150 min solar time after a statis- 
tically stationary state of the magneto-convective system 
had evolved. For the detailed analysis we then considered 
10 equidistant snapshots covering about 20 min in order to 
suppress the effects of the 5-min ute osci l lation s and to re- 
duce the statistical fluctuations. IVoglerl (|2005l ) used these 
simulations to study the centre-to-limb variation of the 
bolometric intensity in dependence of the average vertical 
magnetic field. 



where 



R — /J R4 



(2) 



are the respective numbers of data points. 

The reduced x 2 value (i.e. the x 2 value divided by 
the number of degrees of freedom) between the mean of 
dataset I and dataset II is 0.42, suggesting that they are 
equivalent within the scatter. With similar results ob- 
tained for the other wavelength bands, this justifies the 
use of the same PSF for both datasets. 

For the comparison with simulations (described in the 
next section), an extract (or subimage) of an observation 
in this work is chosen such that it has the same physical 
size as one simulation snapshot, which corresponds to ap- 
proximately 7.5 arcsecx 8.3 arcsec under a limb angle of 
/j, = 0.9. Such an individual segment can then be directly 
compared with one simulation in terms of visual appear- 
ance, histograms, and contrast, as is done in Sect. 4. 

3. Calculations 

MHD simulations: The 3-D MHD simulations have 
been run using the MURaM code (jVogler et alj2005l) with 
non-grey radiative transfer using 4 opacity bins and 24 ray 
direct ions for determinin g the radiative heating/cooling 
rates ( Voeler et al. 20041) . On the Sun, the computational 
box corresponds horizontally to a square of 6 Mm x 6 
Mm and has a vertical extension of 1.4 Mm, covering the 
range between about 500 km above and 900 km below the 
visible solar surface that is located at an average height 
of optical depth unity. The cell size of the numerical grid 
is 20.8 km in the horizontal directions and 14 km in the 
vertical. We used 3 simulation runs with initially homo- 
geneous and unipolar vertical magnetic fields of 0G, 50G, 
and 200G, respectively. Since we wish to study the effect of 
the surface magnetic field on the radiative output at the 
surface, we do not prescribe the total radiative surface 
flux, but rather assume a fixed constant entropy density 
of the fluid enter ing the simul ation box through the open 
lower boundary (jVoglerl 120051 ) . The value of the entropy 
density has been chosen such that the resulting mean ra- 
diative energy flux in the non-magnetic run matches the 
observed solar value to within 1%. In all three simulations, 
the total radiative output at the top of the simulation box 
fluctuates by about ±1% around its mean value in the 



Intensities: We calculate the emergent intensities along 
parallel rays, each passing through a different point of the 
surface grid of the MHD simulations. Since we are inter- 
ested in a broad spectral c overage, we use Kurucz' spectral 
synthesis code ATLAS9 (jKuruczl [1993) (which includes 
diatomic molecules) with opacity distribution functions 
(ODFs) on each of those simulated model atmospheres at 
disk centre and towards the limb. ATLAS9 (as rewritten 
by J. B. Lester) and the ODFs were obtained through 
CCP7 (Collaborative Computing Project No. 70). The 
program solves the radiative transfer equations under the 
assumption of local thermodynamic equilibrium (LTE). 
The simplification of LTE and the use of ODFs allows us 
to calculate the whole spectrum from 160 nm to 160 000 
nm with a spectral resolution of better than 200 in the 
visible. However, the HINODE/SOT/BFI filters are nar- 
rower than the ATLAS9 resolution bins, and, in addition, 
the spectra are not calculated at the wavelength centre for 
the filters which leads to a comparison at slightly different 
wavelengths. Furthermore, the CN and GB regions com- 
prise many spectral lines that are only approximated by 
coarsely resolved ODFs in the ATLAS9 calculation. This 
may influence the comparison of the simulated intensities 
with the observations for the CN and GB wavelengths. 

For an inclined view away from disk centre, we allow 
the line of sight to cut through the simulated atmosphere 
at an angle 8 to the surface normal (0 is the heliocentric 
angle). For our simulation box, this means we have to 
interpolate the atmosphere parameters along rays with the 
respective heliocentric angle. 

For the comparison with observations, we degraded 
the simulated images to take into account optical effects: 
The point-spread fun ctions for the SOT as calculated by 
Mathew et al. ( 20091 ) were used for both dataset I and 
dataset II, as discussed in Sect. 2. 

4. Comparisons of observations with simulations 

Here, we compare the resulting intensities from simula- 
tions with observations to validate the simulated results. 
For dataset I at /x = 0.9 with its finer binning, we first in- 
spect the observed and simulated intensity images visually 
for a qualitative comparison. In a next step we compare 
dataset II with the simulations (at all limb angles) in terms 
of rms contrasts and intensity distributions. 



http:/ /ccp7. dur.ac.uk/ 



N. Afram et al.: Intensity contrast from MHD simulations and from HINODE observations 



5 




Fig. 5. Histograms of the intensity distribution in convolved simulations and in observations for the 5 Hinode/SOT 
wavelengths, from top to bottom: 388 nm, 430 nm, 450 nm, 555 nm, 668 nm. From left to right: Limb angles [i = 0.9, 
\i = 0.6, \i = 0.3. The histograms for the simulations are shown for 3 different magnetic fluxes: 0G (dotted), 50G 
(dashed), 200G (long dashes). To illustrate the spread in the observations, we overplot histograms of extracts of the 
whole observed image (grey lines). The black solid line indicates the mean over all observed extracts. Note the different 
scales in the horizontal axes. 



4.1. Simulations vs. dataset I: Visual comparison 

In this section, we compare observations with original and 
convolved simulations at a limb angle of fi = 0.9. In Fig. [3] 
original and convolved MURaM simulation snapshots for 
50G are shown together with the observed image in the 
respective wavelength filter. The qualitative agreement is 
good, the granules are of comparable size and the visually 



perceived contrasts appear to be similar. The occurrence 
of brighter and darker features is similar and also changes 
correspondingly for the different wavelength filters. The 
brightest features can be seen in the CN wavelength. The 
bright magnetic features that appear for the 50G (as seen 
in Fig. [3]) and 200G simulations are not seen for nonmag- 
netic simulations. 
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4.2. Simulations vs. dataset II: rms contrasts and 
intensity distributions at different limb angles 

In Fig. @] we demonstrate the effect of convolving with 
the PSF in a more quantitative way. We plot, as an ex- 
ample, the histograms for the green continuum intensity 
distribution in the original simulation at the limb angle 
(j, = 0.9 for OG, 50G, and 200G. We overplot (in grey) the 
corresponding histogram for the convolved simulation as 
shown in Figure [5] As expected, the convolved histograms 
are significantly narrower. For most wavelengths we find 
a distribution with a double peak in the unconvolved 
simulations, which is most prominent in the 50G simu- 
lation, and represents the dark intergranular components 
(the more prominent peak at 1/(1) <1) and the bright 
gra nular components (the weaker peak at 1/( 1) >1) (see 
e.g. IWedemever-Bohm fc van der Voortl ( 20091 ) and refer- 
ences therein). After the convolution, the double-peak is 
smeared out, but important differences between the non- 
magnetic, 50G, and 200G simulations remain, in particu- 
lar for the visible broadband distributions (BC, GC, and 
RC). 

Figure [5] shows the histograms of the intensity for (all) 
10 averaged simulation snapshots and for average mag- 
netic fields of OG, 50G, and 200G. From left to right we 
plot histograms for the convolved simulations at the three 
limb angles fi — 0.9, \i — 0.6, fi — 0.3 (i.e. from disk centre 
towards the limb) and for all 5 Hinode/SOT wavelengths 
(from top to bottom: CN: 388 ran, GB: 430 nm, BC: 450 
nm, GC: 555 nm, RC: 668 nm). We also show the averaged 
histograms for the Hinode observations (dataset II, solid 
black line) . We illustrate the spread in the observations by 
overplotting histograms of extracts of the observed image 
(grey lines). For ease of comparison, the simulations have 
been renormalised to cover the same number of points as 
the observations. 

At the limb angle of fi = 0.9, the agreement between 
data and simulations is reasonable for all wavelengths. 
Closer to the limb, at fi = 0.6 and \i = 0.3, the histograms 
for the simulations in the CN and GB wavelengths do not 
agree with the measurements. For the three continuum 
filters, however, the agreement is satisfactory, as the his- 
tograms for the simulations all fall within the scatter of 
the data. None of the simulations exactly matches the ob- 
served histograms, which leads to high x 2 values. A quan- 
titative comparison and the \ 2 analysis are discussed in 
the next section. 

Intensity variations between bright and dark features 
determine the root-mean-square (rms) contrast, which 
is the standard deviation divided by the mean value. 
Figure [6] shows the rms contrasts for the convolved 0G, 
50G, and 200G simulations in comparison to the result- 
ing rms contrast from the observations at the limb angles 
fi = 0.2 — 1.0. For the simulations, these rms contrasts are 
the averages over rms contrasts of all simulation runs. For 
the observations, we calculate the rms contrasts for each of 
the extracts of the observed images, i.e. extracts of similar 
size as one simulation snapshot, and then average them. 
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Fig. 6. Rms contrasts for the BC, GC, and RC 
Hinodc/SOT wavelengths at the limb angles /i = 0.2 — 1.0 
for the mean over the 0G simulation snapshots (dotted 
lines), the 50G simulation snapshots (dashed line), the 
200G simulation snapshots (long dashes) and the Hinode 
dataset II observations (diamonds). We show the spread of 
the 10 snapshots of the 50G simulation and overplot the la 
error of the Hinode dataset II observations. Overplotted 
(black crosses) are the Hinode dataset I results at (i = 0.9, 
which are almost indistinguishable from the dataset II re- 
sults. 



We consider the change in the rms contrast as one 
moves from weaker to higher magnetic fluxes, as is in- 
dicated by the dotted, dashed, and long-dashed lines in 
Fig.|5]for the 0G, 50G, and 200G simulations, respectively. 
It can be seen that the 200G rms contrast is not a simple 
extrapolation of the 0G and the 50G contrasts. The 0G 
and 50G rms contrasts look very similar, although the 50G 
rms contrast falls off a bit more steeply. The 200G rms con- 
trast shows a much flatter profile with a lower contrast at 
disk centre and a higher contrast at the limb compared to 
the weaker fields. The rms contrasts suggest that the 50G 
simulation agrees slightly better with the Hinode observa- 
tions than the 0G simulation for all wavelengths. Although 
the 200G simulation fits the observations better near disk 
centre (as can also be seen from the \ 2 statistic in the next 
section), its overall behaviour does not represent the one 
seen in the data. Given that the observations refer to quiet 
Sun, this result is not surprising. The agreement between 
the contrast as measured from the observations and the 
simulation for 50G is good in the visible continuum bands 
at 450.4 nm (BC), 555.0 nm (GC), and 668.4 nm (RC). 
The contrasts derived for 388.3 nm (CN) and for 430.5 
nm (GB) are generally too large. For the 50G simulation, 
they can exceed the measured contrasts by up to a factor 
of 1.5 and 1.3 for CN and GB, respectively, depending on 
limb angle. 
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Fig. 7. x 2 values for histograms of convolved simulations as shown in Fig. [5] but for all the limb angles fi = 0.2 — 1.0. 
The x 2 values for OG, 50G, and 200G are represented by dotted, dashed, and long-dashed lines, respectively. The A\ 2 
for a confidence level of 99% is 75 for 50 bins. A A% 2 of 75 is indicated by the solid line on the right of each plot. The 



grey line at x = 49 indicates the level where xi = !• 




4.3. Statistics for simulations 

Although the differences between the simulations with dif- 
ferent average magnetic fluxes are less distinct after the 
convolution, they can still be seen in the histograms in 
Fig. [5] To find the simulations that agree best with the ob- 
servations, we calculate the value of x 2 according to Eq. 1, 
where Ri and Si are the number of pixels in intensity bin i 
for the observations and simulations, respectively. The x 2 
values are shown in Fig. [7] for limb angles of/i = 0.2 — 1.0 
and for the five Hinode wavelength filters. As each his- 
togram has 50 bins, i.e., 49 degrees of freedom (d.o.f), we 
would expect a good fit to have a value of x 2 — 49. This 
corresponds to a value of approximately one for the re- 
duced x 2 given by Xi = X 2 /d.o.f. The grey line in Fig. [7] 
indicates the level where Xr = !• 

It turns out that, except at disk centre, none of the 
simulations fit the CN or G-band observations well. This 
is not too surprising as we are currently assuming LTE in 
our spectral synthesis calculations, and are, furthermore, 
not carrying out a full line-by-line synthesis, but use ODFs 
to approximate the opacity in relatively coarse wavelength 
bins (see Sect. 3). For the blue, green and red continuum 
bands, the simulations agree reasonably well, though it 



appears that the 200G simulation fits best for positions 
closer to disk centre, while the OG and 50G simulations 
show better agreement closer to the limb. 

We can use contours of equal x 2 to examine whether 
the convolved intensity distributions resulting from differ- 
ent magnetic flux levels are distinguishable, i.e., appear to 
be drawn from a different underlying population. To this 
aim, we calculate A% 2 = X 2 — xL where Xm denotes the 



minimum x 2 ■ I n our case of 50 bins, A% 2 values in excess 
of 75 sugge st a different distr ibution at a confidence level 
of 99% (see lPress et~al]|l992h . The 99% confidence levels 
are indicated by the vertical bars on Fig. [7] They con- 
firm that the intensities calculated for the CN filter are 
not a good representation of the observations, except near 
disk centre, where the convolved non-magnetic, 50G and 
200G distributions are indistinguishable at the 99% confi- 
dence level. In all other bands, the non-magnetic and the 
50G convolved distributions are largely indistinguishable, 
while the 200G distributions appear to show significantly 
different behaviour, in particular towards the limb (see 
also the right-hand column in Fig. [5]). While the 50G sim- 
ulation does not represent the best fit at most limb angles, 
it largely lies within Ax 2 = 75. 
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5. Intensity contrasts in simulations 

The centre-to-limb variation of the emergent intensity 
varies with wavelength. In the UV, magnetic features are 
brighter than at other wavelengths. These differences be- 
come even stronger near the limb and are therefore crucial 
for spectral solar irradiance variations. It is generally not 
possible to reproduce the observed limb behaviour with 
1-D model calculations, as they tend to overestimate the 
limb contrast at extreme limb angles. 

In Figure [5] we show the simulated contrasts (i.e. the 
mean intensity of the 10 averaged magnetic 3-D MHD 
simulation snapshots compared to the mean quiet Sun 
(= 0G) intensity) with average vertical magnetic field of 
50G and 200G as a function of heliocentric position in 
the five Hinode wavelength filters. As expected, the con- 
trast increases towards the limb and is higher for shorter 
wavelengt hs. The resul t is qualitatively consistent with 
figure 4 in IVogl er (2005) if one assumes that the bolomct- 
ric intensity corresponds roughly to the maximum of the 
Planck curve. Contrary to 1-D models, the flattening and 
turnover (e.g. for 50G at 388 nm and 430 nm) are in better 
agreement with measurements of facular contrasts than 
1-D calculations (|Topka et al.l Il997t lUnruh et all 1 19991 
Ortiz et alj|2002h . 



6. Conclusions and Outlook 

In the present paper, we have compared simulated and 
observed intensity distributions in the quiet Sun at a range 
of limb angles and derived mean contrast for simulations 
with an average magnetic field of 0G, 50G, and 200G. 

Comparing histograms and rms contrasts of the ob- 
servations and simulations, we find that no single simula- 
tion reproduces the data for all limb distances and colour 
bands. The 50G simulation intensity histograms and rms 
contrasts agree best with the observations near the limb 
and show an overall similar behaviour. The 200G simula- 
tion, however, reproduce the observations better near disk 
centre, but do not represent the overall behaviour of the 
observed rms contrasts. 

There may be various reasons for the discrepancy be- 
tween the data and the simulations. One of these could be 
the use of LTE in ATLAS9 along with the fact that we 
are not calculating spectral lines (which mainly affects the 
results for CN and GB). Also, in our dataset we are observ- 
ing different areas on the Sun at various limb angles with 
potentially differing activity levels. In additon, by aver- 
aging over a large area, the Hinode histograms represent 
a de facto mixture of features/magnetic fluxes, that we 
would expect to be reproduced by a combination of the 
simulations with 0G/50G/200G average magnetic fields. 
Furthermore, the height of the model atmosphere in our 
simulation box might be insufficient. As the radiation is 
produced at increasingly higher layers in the atmosphere 
as one moves to the solar limb and to shorter wavelengths, 
the computational box might not be sufficiently high to 
encompass enough of the region where the radiation origi- 



nates. This problem becomes particularly acute for the CN 
and GB data, due to the presence of strong lines formed 
particularly high. 

The approach presented in this study should allow us 
to calculate the facular and network contrast not solely as 
a function of limb angle, but also as a function of magnetic 
flux. For this purpose, we will calculate magnetograms 
from MURaM simulations to identify the previously deter- 
mined emergent intensities with corresponding magnetic 
fluxes. Depending on further tests, these facular and net- 
work intensities could then be incorporated into irradiance 
reconstruction models (such as SATIRE; see e.g. Krivova 
et al. 2003) which would enable us to help improve these 
models, as well as identify the cause of irradiance varia- 
tions even more reliably. Of particular interest will also be 
a detailed comparison of the output of MHD simulations 
with t he high-resolution, low-stray light Sunr ise observa- 
tions ( Barthol et al. 20101 : Solanki et al. 2010l). especiall y 
those from the SUFI instrument ( Gandorfer et al.|l2010h . 
which will go beyond j ust the comparison of rm s values at 
solar disk centre as in iHirzberger et ah ( 2010l ). 
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